Astronomy & Astrophysics manuscript no. off ' accepted 


©ESO 2010 


February 23, 2010 





Ultraluminous X-ray sources out to z~ 0.3 in the COSMOS field. 

V. Mainieri 1 , C. Vignali 2 - 3 , A. Merloni 4 - 5 , F. Civano 6 , S. Puccetti 7 , M. Brusa 5 , R. Gilli 3 , M. Bolzonella 3 , A. Comastri 3 , 
G. Zamorani 3 , M. Aller 9 , M. Carollo 9 , C. Scarlata 9 - 10 , M. Elvis 6 , T.L. Aldcroft 6 , N. Cappelluti 5 , G. Fabbiano 6 , 
A. Finoguenov 5 , F Fiore 8 , A. Fruscione 6 , A.M. Koekemoer 11 , T. Contini 12 , J.-P. Kneib 13 , O. Le Fevre 13 , S. Lilly 9 , 
A. Renzini 14 , M. Scodeggio 15 , S. Bardelli 3 , A. Bongiorno 5 , K. Caputi 9 , G. Coppa 3 , O. Cucciati 16 , S. de la Torre 13 , 
L. de Ravel 13 , P. Franzetti 15 , B. Garilli 15 , A. Iovino 16 , P. Kampczyk 9 , C. Knobel 9 , K. Kovac 9 , F. Lamareille 12 , J.-F. Le 
Borgne 12 , V. Le Brun 13 , C. Maier 9 , M. Mignoli 3 , R. Pello 12 , Y. Peng 9 , E. Perez Montero 12 , E. Ricciardelli 17 , 
J. D. Silverman 9 , M. Tanaka 1 , L. Tasca 13 , L. Tresse 13 , D. Vergani 3 , E. Zucca 3 , P. Capak 19 , O. Ilbert 13 , C. Impey 18 , 

M. Salvato 19 , N. Scoville 19 , Y. Taniguchi 20 , J. Trump 18 



ESO, Karl-Schwarschild-Strasse 2, D-85748 Garching bei Miinchen, Germany 

Dipartimento di Astronomia, Universita degli Studi di Bologna, Via Ranzani 1, L40127 Bologna, Italy 
INAF-Osservatorio Astronomico di Bologna, Via Ranzani 1, 1-40127 Bologna, Italy 
Excellence Cluster Universe, TUM, Boltzmannstr. 2, D-85748 Garching bei Miinchen, Germany 
Max-Planck-Institute fur Extraterrestrische Physik, Postfach 1312, 85741, Garching bei Miinchen, Germany 
Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 USA 
ASI Science Data Center, via Galileo Galilei, 00044 Frascati Italy 

INAF-Osservatorio astronomico di Roma, Via Frascati 33, 00040 Monteporzio Catone, Italy 

Institute of Astronomy, Swiss Federal Institute of Technology (ETH Honggerberg), CH-8093, Zurich, Switzerland 

Spitzer Science Center, Pasadena, CA, 91 125 

Space Telescope Science Institute, Baltimore, Maryland 21218, USA 

Laboratoire d'Astrophysique de Toulouse-Tarbes, Universite de Toulouse, CNRS, 14 avenue Edouard Belin, F-31400 Toulouse, 
France 

Laboratoire d'Astrophysique de Marseille, Marseille, France 

INAF-Osservatorio astronomico di Padova, Vicolo DeH'Osservatorio 5, 1-35122 Padova, Italy 

DSfAF - IASF Milano, Milan, Italy 

INAF Osservatorio Astronomico di Brera, Milan, Italy 

Dipartimento di Astronomia, Universita di Padova, Padova, Italy 

Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721 

California Institute of Technology, MC 105-24, 1200 East California Boulevard, Pasadena, CA 91 125 UA 

Research Center for Space and Cosmic Evolution, Ehime University, Bunkyo-cho 2-5, Matsuyama 790-8577, Japan 



Received accepted ... 

ABSTRACT 

Context. Using Chandra observations we have identified a sample of seven off-nuclear X-ray sources, in the redshift range z=0.072- 
0.283, located within optically bright galaxies in the COSMOS Survey. All of them, if associated with their closest bright galaxy, 
would have L[0.5 - 7 keV] > 10 39 erg s _1 and therefore can be classified as ultraluminous X-ray sources (ULXs). 
Aims. Using the multi-wavelength coverage available in the COSMOS field, we study the properties of the host galaxies of these 
ULXs. In detail, we derived their star formation rate from Ha measurements and their stellar masses using SED fitting techniques 
with the aim to compute the probability to have an off-nuclear source based on the host galaxy properties. We divide the host galaxies 
in different morphological classes using the available ACS/HST imaging. 

Methods. We selected off-nuclear candidates with the following criteria: 1) the distance between the X-ray and the optical centroid 
has to be larger than 0.9", larger than 1.8 times the radius of the Chandra positional error circle and smaller than the Petrosian radius 
of the host galaxy; 2) the optical counterpart is a bright galaxy (Rab <22); 3) the redshift of the counterpart is lower than z= 0.3; 
4) the source has been observed in at least one Chandra pointing at an off-axis angle smaller than 5'; 5) the X-ray positional error is 
smaller than 0.8". We verified each candidate super-imposing the X-ray contours on the optical/IR images. We expect less than one 
misidentified AGN due to astrometric errors and on average 1.3 serendipitous background source matches. 

Results. We find that our ULXs candidates are located in regions of the SFR versus M* plane where one or more off-nuclear detectable 
sources are expected. From a morphological analysis of the ACS imaging and the use of rest-frame colours, we find that our ULXs 
are hosted both in late and early type galaxies. Finally, we find that the fraction of galaxies hosting a ULX ranges from « 0.5% to 
ss 0.2% going from Lo.5-2 kev = 3 x 10 39 erg s _1 to Lo.5-2 w = 2x 10 40 erg s -1 . 

Key words. Xrays: galaxies - X-rays: binaries - X-rays: general - Surveys 



1. Introduction 

An intriguing class of X-ray objects are the so called ultra- 
luminous X-ray sources (ULXs). Here an ULX is defined as 



an X-ray source in an extra-nuclear region of a galaxy with 
an observed luminosity in excess of 10 39 erg s in the 0.5- 
7 keV band. Such X-ray luminosities are higher than expected 
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Fig. 1. The distribution of the difference between the detected 
X-ray positions and the input positions in units of the X-ray po- 
sitional error. 



for spherical Eddington-limited accretion onto a ~ 1OM black 
hole. ULXs were know n already from studies with Einstein , 
ROSAT, and ASCA (e^ . lFabbianoll989HColbert & Ptakl2002t 
Makishim a et al. | |2000), but it was after the advent of Chandra 
with its combination of high angular resolution and moderate 
spectral resolution that has been possible to make significant 
progr esses in their study (e.g. lRoberts e t al. 20041: [Swartz et al.l 
2004). There is a wide debate in the literature on the nature of 
these sources. ULXs may be powered by accretion onto stellar- 
mass black holes assum ing that there is relativistic be aming 
(e.g. Kording et al."l 120021) . or radiative anisotropy (e.g. King 
2002 ), or they may b e associated with super-Eddington discs 
(e.g. iBegelman 1 120021) . It has been also suggested that ULXs 
represent a ne w class of intermediate-mass (10 2 - 10 5 M Q ) 
black holes (e.g. lColbert & Mushotzkvlll999blMiller & Colbert! 
|2004|) . These intermediate-mass black holes may be fed by 
Roche lobe overflow from a tidal c aptured stellar comp anion 
that is not destroyed by tidal heating dHopman et al. 1 12004). Off- 
nuclear AGN activity could be also a signature of a recoiling 
massive black hole: a massive black hole binary coalesces and 
gives origin to gravitational waves which can give a kick to the 
center of mass of the system. If the recoiling black hole retains 
the inner parts of its accretion disk, we could see its luminous 
phase as an off-nuclear AGN (see IVolonteri & Madaul 120081 
and references therein). Finally, ULXs could also be the high- 
luminosity extension of supernovae (e.g. lSwartz et al. 1 12004). 
Many of the previous studies based on Chandra data are focused 
on local galaxies, where the Chandra angular resolution allows 
to detect several off-nuclear sources in one single galaxy. In this 
paper, we select a sample of ULXs from the Chandra survey 
in the COSMOS field. We have here the advantage to combine 
deep X-ray observation with a wealth of multiwavelength an- 
cillary data that we will use to put constraints on the nature of 
these sources and on the properties of their host galaxies. The 
redshift range that we are covering is up to z — 0.3. A study on 
off-nuclear sources i n a similar redshift range was performed by 
iLehmer et al~1 d2006l) on the Chandra Deep Fields (CDFs). 
We quote in this paper magnitudes in the AB system and we as- 



Fig. 2. X-ray to optical offset s in arcsec fo r X-ray sources with 
a secure identification (ICivano et al. I l2010l) and with X-ray po- 
sitional error smaller than 0.8". The circle of 0.9"radius encom- 
passes 95% of the X-ray sources. 

sume a cosmology with Ho = 70 km s _1 Mpc -1 , Qm = 0.3 and 
Q A = 0.7. 

2. Sample selection 

We have selected off-nuclear X-ray candidates from the Chandra 
COSMOS Survey (C-COSMOS), which is a recently completed 
1.8 Ms Chandra program to image the central 0.9 deg 2 of the 
COSMOS field with an effective exposure ranging from -160 
ksec to ~ 80 ksec goin g from the center to the borders of the 
field (Elv is et aTl 12009). The limiting source detection depths 
are 1.9 x 10~ 16 erg cirT 2 s" 1 in the [0.5-2 keV] band, 7.3 x 10~ 16 
erg cirT 2 s _1 in the [2-10keV] band, and 5.7x 10~ 16 erg cirT 2 s _1 
in the [0.5-10 keV] band. We used a point source catalogue in- 
cluding 1761 objects detected in at least one band (0.5-2, 2-7 and 
0.5-7 keV) with a maximum likelihood ratio larger than detml= 
10.8, corresponding to a probability of ~ 2 x 10~ s that a catalo g 
source is instead a background fluctuation dPuccetti et al. 12 009). 
The optical and infrared i dentifications of alm ost all (99.7%) of 
the sources are reported in lCivano et al. I d201 0fT. 

As a first step to select off-nuclear X-ray sources, we veri- 
fied the X-ray position accuracy that we have in the C-COSMOS 
ob servations foll o wing the procedure presented in Sec. 4.3 
of iPuccetti et all (§009). A set of 49 Chandra ACIS-I point- 
ings has been simulated with the MARX0 simulator, adopt- 
ing the same exposure times, aim points, and roll-angles as 
the real C-COSMOS p ointings. The detection code PWDetect 
(iDamiani et al. L 1 19971) was applied to the simulated data. We 
then compared the output of the detection algorithm with the 
input catalogue of the simulation. In Fig. [1] we show the dis- 
tribution of the difference between the detection algorithm po- 
sitions and the input positions in units of the X-ray positional 
error. The latest has been estimated as the ratio of the PSF at 

1 The ULX cand idates presented in this paper are flagged as 'off- 
nuclear' sources in ICivano etal.U2010l) . 

2 http://space.mit.edu/CXC/MARX. 
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Fig. 3. Cutouts in the HST/ACS F814W band dKoekemoer et al. II2007I) of the seven X-ray off-nuclear sour ces in the C- COSMOS 
field. The red cross indicates the position of the X-ray centroid and the red circle the X-ray positional error dElvis et al, 1 l2009h . We 
provide for each object: the Chandra ID (top-left), the redshift (top-right), the logarithm of the X-ray luminosity in the [0.5-7] keV 
band (bottom-left), the maximum likelihood ratio for the X-ray detection (bottom-right), the morphological classification of the 
host galaxy (bottom-middle; see Section [3TTb . The images have different sizes for display purposes; the vertical bar in each cutout 
corresponds to 2" . On the right of each ACS cutout, there is the corresponding Chandra [0.5-7 keV] image. 



the position of the source and the square root of the net, back- 
gr ound subtracted, sourc e counts. In comparison with Fig. 10 
of iPuccetti et al. I ([2009), we have restricted the analysis only 
to sources that have been detected at least in one image at an 
off-axis angle smaller than 5' to take advantage of an excel- 
lent PSF. From the distribution in Fig.[TJ we find that 94% of the 
sources have offsets below 1.8 times the positional error. We will 
adopt this value as a threshold to select off-nuclear candidates 
and therefore we expect that up to 6% of our sample is contam- 
inated by nuclear X-ray sources with large astrometric errors. 
We will shortly come back to this issue. Another possible source 
of spurious off-nuclear objects could be a poor astro metric ac- 
curacy of the X-ray images. According to Fig. 6 of lElvis et al. I 
(2009), 95% of the Chandra sources have an absolute astromet- 



ric accuracy better than 1.4". For our study we aim at even better 
astrometric accuracy, therefore we have considered only the X- 
ray sources with an X-ray positional error smaller that 0.8". We 
show the comparison between X-ray coordinates and optical co- 
ordinates for sources with a secure identification in Fig. [2] 95% 
of the X-ray sources have an absolute astrometric accuracy bet- 
ter than 0.9". 

Summarizing, the off-nuclear candidates were selected using 
the following criteria: 

a) The distance between the X-ray centroid and the optical cen- 
troid has to be larger than 1 .8 times the radius of the Chandra 
positional error circle at that position. 

b) The X-ray positional error is smaller than 0.8". 
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Fig. 4. Aperture photometry for each off-nuclear candidate. The vertical dot-dashed line indicates the distance between the X-ray 
position and the centroid of the host galaxy; the v ertical dashed line is t he Petrosian radius of the host galaxy; the horizontal line 
corresponds to the counts estimated by EMLdetect dPuccetti et al. LI2009I) . For objects XID= 1151, 1388, and 1870, the filled circles 
represent the photometry on an area including 90% of the PSF obtained applying an aperture correction factor to the photometry 
measured on an aperture of 1 "radius (see discussion in Sec.|2]i. 



c) The source has been observed in at least one Chandra point- 
ing at an off-axis angle smaller than 5'. 

d) The optical counterpart is a bright galaxy (Rab <22). 

e) The redshift of the host galaxy is less than z=0.3. The pro- 
jected linear distance corresponding to an average Chandra 
positional error is ~ 4 kpc at z=0.3. This means we will con- 
sider only off-nuclear candidates that are more than ~ 7 kpc 
away from the center of the galaxy at z=0.3. At larger red- 
shifts we would be able to select only off-nuclear candidates 
that are at larger distances (> 7 kpc) from the host galaxy 
center, where the n umber of observed off-nuclear sources 
seems to decrease dSwartz et al. 1 12004b and we would be 



more affected by the contamination of background objects. 
Therefore we limit our sample to z< 0.3. 
f) The distance between the X-ray centroid and the optical cen- 
troid is larger tha n 0.9"and smaller than the Petrosian radius 
dPetrosian II 19761 RrEI) of the galaxy, which we use as a mea- 
sure of the galaxy's extension. 



3 Rp is defined as the radius at which the ratio (rp) of the local sur- 
face brightness at that radius and the mean surface brightness within 
that radius equals some specified value rpjj m . For a surface brightness 
distribution described by a de Vaucouleurs or an exponential profile, 
a value rpjj„, = 0.2 is reached at Rp ~ 1.8 R1/2 and Rp ~ 2.2 R1/2, 
respectively (R1/2 is the half-light radius of the galaxy, see Fig. 17 of 
IScarlata et all2007h . 
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If we consider the selection criteria b), c), d), and e), only 
sixteen sources from the C-COSMOS catalogue satisfy all of 
them. Based on our previous discussion of Fig. Q] we expect up 
to 6% of our sample to be due to nuclear X-ray sources with large 
astrometric errors, therefore we can conclude that our sample of 
ULXs contains less than one misidentified AGN. 

For all the candidates provided by these selection criteria we 
have verified that no other counterpart closer to the Chandra po- 
sition was present in any band from the u* {A center = 374.3 nm) 
filter to 24 micron. After this one-by-one check, we were left 
with seven off-nuclear source candidates. Cutouts of these ob- 
j ects, obtained from the COSMOS HST/ACS F814W imaging 
dKoekemoer et al. 112007}) . are shown in Fig. [3] together with the 
corresponding Chandra [0.5-7 keV] image. 

Each one of our off-nuclear sources has an estimate of 
the X- ray flux in the [0.5-7] keV band reported in lElvis et aTl 
(l2009h . These fluxes are derived from the counts estimated by 
EMLd etecffl co rrecte d to an area including 90% of the PSF 
dPuccetti et al. L l2009h . In some cases such an area is large 
enough to include the whole host galaxy and therefore the X- 
ray flux could be the total integrated flux of the host galaxy it- 
self. This would include the contribution from the population of 
X-ray binaries in the host, emission from diffuse gas and a pos- 
sible weak central AGN. In order to estimate these possible con- 
taminations on the measured X-ray fluxes, we have performed 
aperture photometry for each off-nuclear source. The radii of the 
apertures have been chosen with increasing size from a mini- 
mum of 1" up to include the whole galaxy. In Fig. |U we plot 
the net counts in the [0.5-7] keV band as a function of the aper- 
ture radius. For four of our sources (XID= 2418, 3441, 11100, 
11938) the counts measured at different apertures are constant 
within the uncertainties. Therefore, we assume that the contri- 
bution of the host galaxy is not significant compared with the 
uncertainties on the measure. For the remaining three sources 
(XID= 1151, 1388, 1870) the counts rise with the aperture ra- 
dius and there may be a significant contamination due to the 
integrated flux of the whole galaxy. In order to minimize this 
contamination, we have considered the measured counts in the 
smaller aperture (1")- We have then used the known PSF shape 
at the position of the source to estimate the expected fraction be- 
tween the counts measured in an aperture of l"and the ones over 
an area corresponding to 90% of the PSF. We have then used this 
ratio to convert our measured counts on the 1 "aperture into the 
expected ones on a 90% PSF area. These corrected counts are 
indicated with a filled circle in the plots of Fig. |H and we have 
used them in order to estimate the X-ray fluxes. 
Full band 0.5-7 keV fluxes and errors were computed converting 
counts rates to fluxes using the formula: Flux=B ra , c /(CF*10 n ), 
where B rale is the count rate estimated as described above, and 
CF is the energy conversion factor. This conversion factor varies 
with the energy band and the spectral index F assumed for 
the power-law spectrum. We have used the corre ction factor 
CF= .89 counts erg 1 cm 2 reported in Tab. 4 of lElvis et al. I 
(2009) obtained for the 0.5-7 keV band and T = 1.7. We decided 
fo r this average value of the spectral index following the study 
of ISwartz et aU d2004l) that has found a mean power-law index 
of r = 1 .74 + 0.03 for a sample of 154 ULX candidates observed 
with Chandra. We finally report on Tab. Q] the 0.5-7 keV lumi- 
nosities and errors for the seven off-nuclear sources. All sources 
have luminosities well in excess of 10 39 erg s _1 in the 0.5-7 keV 
band (the lowest X-ray luminosity in this band is * 9 x 10 39 erg 
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Fig. 5. X-ray luminosity in the 0.5-7 keV band vs. redshift of 
the seven off-nuclear sources. The different symbols correspond 
to the host galaxy classification based on morphology and rest- 
frame colours: circles are ETGs, while stars are LTGs (see Sec. 
13. II). S quares are the off nuclear sources from iLehmer et al. I 
(2006 |); crosses are the collection of local off-nuclear sources 
by Liu & Mirabel (2005). The dashed line corresponds to the 
flux limit in the deepest region of the C-COSMOS survey: 
SumtO.5-7] =4.7x 10- 16 cgs. 



s -1 ) and are therefore classified as ULX sources, using either 
spectroscopic or photometric redshifts. 

We have secure spectroscopic redshifts for four host galax- 
ies fr om zCOSMOS VIMOS observations at VLT dLillv et al. I 
2007, 2009). For the remaining three objects we have used 
the extremely accurate photometric redshifts availa ble in the 
COSMOS field dllbert et al. l2009HSalvato et al. 120091) based on 
30 broad, intermediate, and narrow bands from the UV to the 
mid-IR. We show in Fig. [5] the X-ray luminosity in the [0.5-7] 
keV band versus redshift of the seven ULXs. The X-ray lumi- 
nosities were computed according to the formula: 



L x = 4nd 2 L fx(l + if' 2 



(D 



http://xmm.esac. esa.int/sas/8.0.0/emldetect 



where dh is the luminosity distance, fx is the X-ray flux 
in the [0.5-7] keV band, and F is the X-ray photon index. We 
assumed r = 1.7 (see discussion in this section). Different 
symbols correspond to the morphological classes of the host 
galaxies (see S ec. 13. U . Squares are the off-nuclear sources from 
Lehmer et al. (2006); crosses are the collection of local ULXs 
by iLiu & Mirabel I (120051) . The dashed line corresponds to the 
flux limit of the C-COSMOS survey, Si im [0.5 - 7] = 4.7 x 10~ 16 
erg cirT 2 s _1 . 

In order to estimate how many background sources we ex- 
pect to contaminate our sample, we applied a random shift be- 
tween 30"and 2' to the C-COSMOS sources and searched for 
chance coincidences with Rab <22 and z< 0.3 galaxies. We re- 
peated this procedure 10,000 times and we found that, on av- 
erage, the chance coincidences are « 1.3. Only for 2% of the 
10, 000 simulations we found more than 3 chance coincidences. 
Summarizing, we expect less than one misidentified AGN due to 
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Fig. 6. Colour-mass diagram: circles and stars are respectively 
ULX host galaxies classified as ETGs and LTGs based on their 
morphology/colours; the dots are galaxies in the C-COSMOS 
area with z< 0.3 and Rab < 22. For the source XID=1151 the 
photometric coverage is limited to few bands and we cannot con- 
strain its stellar mass. 



astrometric errors and on average 1.3 serendipitous background 
source matches. 



vided into three ETGs and four LTGs. We will describe in Sec. 
I3.2l the method used to estimate stellar masses. 

The slight preference for ULXs to be hosted in LTGs could 
be explained by the different shapes of the X-ray luminosity 
function (XLF) for Low-Mass X-ray Binaries (LMXBs) and 
High-Mass X-ray B i naries (HM XBs) derived for local galaxies 
(Grim m et al. Il2003h lGilfanovll2004h : the former has an abrupt 
cut-off at Lx ~ 10 39 erg s _1 , while the latter can be described 
with a power-law with a slope a = 1/6. Since early-type stars 
are the dominant stellar population of LTGs, we expect X-ray 
binaries with O or B type companions, HMXBs, to be common 
in these objects. This translates into a higher chance to detect 
ULXs in LTGs or, in any case, in galaxies with current star for- 
mation activity. 

3.2. Stellar masses and SFR 

Stellar masses (M+) are derived from the stellar population syn- 
thesis model that represents the best fit to the observed photom- 
etry (from the u* band to 4.5 jum) using ax 2 minimization tech- 
nique. The procedure is explained in detail bv lBolzonella et al71 
(120091) . Here we only recall the basic ingredients of the Spectral 
Energy Distribution (SED) fitting procedure: 

- stellar population s ynthe sis models from the libraries of 
iBruzual & Chariot I (12003b : 

- eleven "smooth" star formation histories for each library: 
one constant star formation model plus 10 r-model with e- 
folding time-scales r = 0.1,0.3,0.6, 1,2,3,5, 10, 15,30Gyr; 

- a Chabrier initial mass function; 

- a Calzetti extinction law with < Ay < 3; 

- solar metalicity (Z=Z ). 



3. Host galaxy properties 

3.1. Galaxy classification 

Studies of local samples of ULXs (e.g. lSwartz et al. | |2004) have 
shown that these sources are mainly present in late type galaxies. 
A visual inspection of Fig. [3] suggests that the ULXs at interme- 
diate redshifts that we are studying are hosted in both early and 
late type galaxies (ETGs and LTGs, hereafter). 

To confirm this impression we classified t he host galax- 
ies ba sed on their morphology and colours (e.g. Mi gnoli et al. I 
l2009h . Taking advantage of the COSMOS HST/ACS F814W 
images ( Koekemoer et al. 2007), we h ave used an a c curate 
morphological classification derived by IScarlata et al. I d2007h 
through the Zuric h Estimator of Structural Type (ZEST). 
IScarlata et al. I (l2007b describe in detail the methodology and 
the performances of this method. We only recall here that the 
ZEST classification is based on: a) five non-parametric diagnos- 
tics (asymmetry A, concentration C, Gini coefficient G, 2 nd order 
moment of the brightest 20% of galaxy pixels M20, ellipticity e); 
and b) the exponent n of single Sersic fits to the two-dimensional 
surface brightness distributions. ZEST assigns to each galaxy a 
morphological type (l=early type; 2=disk; 3=irregular) and a 
bulgeness parameter that splits the disk galaxies in four sepa- 
rate bins, from bulge dominated disks (2.0) to pure disk galaxies 
(2.3). For the bulge-dominated galaxies (2.0), we complemented 
the morphological information with their rest-frame colours to 
further subdivide them: if they have red U-B rest-frame colours 
we include them in the ETGs sample (XID= 1151), otherwise 
we classify them as LTGs (XID= 1870). In Fig. [6] we plot the 
colour-mass diagram for our ULX host galaxies: they can be di- 



The star formation rate (SFR) values have been estimated 
using when possible (for three U LXs) the Ha- /16563 line fl ux as 
measured by the routine Platefit (lLamareille et al. 1 12009). with 
a correction for reddening. We used the iKennicutt I (0,998) rela- 
tion between Ha and SFR: SFR(M yr" 1 ) = (7.9 x lO' 42 )L(Ha) 
ergs s _1 . The de-reddened flux of Ha has been computed ac- 
cording to the formula: F der - F obs x 10 c[1+/w] where f(X) = 
3.15854 x lO" 1 - 02109 - 1 - 1 and c = \AlE B -v dSeaton I Il979t 
Mai er et al. 1 12005). If a measure of the W/3 flux was available 
we have estima t ed Eb -v from the Balmer decrement, adopting 
the |Q' Donnell I (1 1 994b Milky Way extinction curve. Otherwise, 
we used the average v alue < Eb-v >~ 0.2 mag derived by 
Moustakas et al. (2006). For the sources without Ha in the spec- 
tral range or for which we have only a photometric redshift, we 
have used the SFR estimate from the SED fitting procedure. 

From the COSMOS catalogue (ICapak et al. I 120071; 
lllbert et al. I g009), we have selected a comparison sample 
of galaxies inside the area covered by Chandra. We imposed 
the same constraints used to select the off-nuclear candidates: 
z<0.3 and Rab < 22. We have also remo ved all the source s that 
are best fitted by stellar SED templates (lllbert et al. 112009b . At 
the end, the comparison sample consists of 2066 galaxies. For 
all of them we have derived stellar masses and SFR values as 
described above. 

We now estimate the probability to have an off-nuclear 
source given a host galaxy with a particular stellar mass and 
SFR. We will consider both LMXBs and HMXBs. 

Fo r LMXBs we used the average XLF derived by Gilf anovl 
(l2004b . This is described by a powerlaw with two breaks, from 
their formula (8): 



V. Mainieri: ULXs in COSMOS 



7 
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Fig. 7. SFR versus stellar masses of the galaxies in the compari- 
son sample (see text). The contours corresponds to region where 
more than 0.1 (red), 1 (green), 5 (cyan), and 10 (yellow) X-ray 
off-nuclear sources per galaxy are expected. The symbols show 
the location in this plane of the host galaxies of the ULXs. For 
the source XID= 1151 the photometric coverage is limited to few 
bands and we cannot constrain its stellar mass. Symbols are the 
same of Fig. [5] The dashed line is the locus where we expect the 
same number of LMXBs and HMXBs with L x > 10 38 erg s _1 . 
Above this line the number of HMXBs is expected to be higher 
than that of LMXBs. 



dN 
dL 3S 



K\ (L3$/Lt,,i) m L38 < Lb t \ 

^2(^38/^6,2) " 2 Lj,,i < L38 <Lb,i 

K3 (L3%/L cut ) a3 Lfc,2 < ^38 < L cut 

L38 > L cu , 



(2) 



where L 3 % 
fined as: 



Lx/10 38 erg/s and normalizations #1,2,3 are de- 



K 2 = 
K 3 



- K2 (Li,,2/L cut ) ai 



We used the best fitting parameter derived by Gilf anov I 
J200l : ai = 1.0, L b ,i = 0.19, a 2 = 1.86, L w = 5.0, a 3 = 4.8. 
The high-luminosity cut-off was fixed at L cut = 500. For the 
average norma lization we used the best fitting value given by 
iGilfanovl (120041) . K x = 440.4 + 25.9 per 10 11 M Q , and we will 
assume a linear relation between the number of X-ray sources 
and stellar mas s as found by the same authors (see Sec. 5 of 
Gilfa nov 1 120041). We note t hat up to L x ~ 2 x 10 39 ergs s" 1 
the XLF of Gilfanov (2004) is consiste nt with later studies (e.g. 
see Fig. 14 of|Humphrev & Buote 2008). Above this luminosity 
we are extrapolating the XLF since no data are currently avail- 
able and therefore the uncertainties are large. For the slope at 
the highest luminosities, we have considered values in the range 
a 3 = [2, 6] and we did not find any significant difference from 
the contours reported in Fig. [7] 



For the HMXBs we used instead the luminosity function de- 
rived by Grimm et al. (2003). In particular, we used the cumu- 
lative form of it, corresponding to their formula (7): 



N(> L 3 s) = 5ASFR(l\ 



-o.6i_2lo-°- 61 



(3) 



where SFR is in units of M yr~' . 

We have then calculated the number of X-ray binaries with 
Lx > 10 38 erg s _1 that we expect in each galaxy integrating the 
XLFs for a given SFR and M*. In Fig. Qwe show the contours 
corresponding to regions where we expect more than 0. 1 (red), 1 
(green), 5 (cyan), 10 (yellow) X-ray sources with Lx > 10 38 erg 
s _1 . In reality these numbers have to be considered upper limits 
because we have not taken into account the limited Chandra spa- 
tial resolution that does not allow to detect off-nucle ar sources 
with small offsets (see Fig. 10 of lLehmer et al. 1 12006). 

From Fig. |7] we find that all our ULX candidates are hosted 
in galaxies for which a large number of X-ray binaries is pre- 
dicted. The dashed line in Fig. [7] is the locus where we expect 
the same number of LMXBs and HMXBs with L x > 10 38 erg 
s _1 . This line clearly divides a region (below the line) where the 
XLF of LMXBs is dominating and therefore the contours are 
mainly defined by the M* values, from a region (above the line) 
where the HMXBs are more numerous and the contours are de- 
termined by the level of the SFR. Our morphological classifica- 
tion is consistent with this picture: ETGs, characterized by lower 
SFR and high stellar masses, are located in the bottom-right part 
of the plot, where the expected number of LMXBs is higher 
than the number of HMXBs. However, we note that there are 
sugge stions in the l iteratu re that no ULX LMXBs may actually 
exist. Ilrwin et aL 1 (12001 had shown that the number of ULXs 
detected in a sample of 28 ellipticals observed with Chandra is 
equal to the number of expected foreground/background objects. 
Additionally, such ULXs are uniforml y distributed and d o not 
follow the optical light of the galaxies. Ilrwin et aT71 ([2004) also 
verified that the same statements can be mad e for the ULXs 
assoc iated to early-type galaxies presented in Colbert & Ptak 
(120021) . 

It would be interesting to repeat the same computation that 
generated Fig.|7]considering only X-ray binaries with Lx > 10 39 
erg s _1 , and therefore to be able to verify the hypothesis that 
ULXs are the high-luminosity tail of normal X-ray binaries. 
Unfortunately, this is not possible due to the poor-knowledge of 
the high luminosity slope of the XLF for LMXBs and HMXBs. 



4. Fraction of galaxies hosting a ULX 

As already pointed out by Ptak & Colbert I d2004), useful con- 
straints on the nature of ULXs can be obtained deriving the frac- 
tion of galaxies that har bor a ULX as a funct ion of the X-ray lu- 
minosity. For example. lKording et al. I (2002) have compared the 
luminosity distribution of X-ray point sources in nearby galax- 
ies with that predicted by X-ray population synthesis models 
to check whether microblazars (microquasars with relativisti- 
cally beamed jets pointing towards the observer) may represent 
an alternative to the intermediate mass black holes scenario for 
ULXs. In order to compute this fraction, we have used the com- 
parison sample selected in Sec. 13.21 We derived for each indi- 
vidual galaxy a 90% upper limit on its X-ray flux in the [0.5-2] 
keV band according to the procedure described in Sec. 6.5 of 
iPuccetti et al. I d2009l) . to which we refer the reader for details. 
The top panel of Fig. [8] shows the number of galaxies for which 
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Fig. 8. Top panel: the number of galaxies for which we could 
detect an off-nuclear source of a given 0.5-2 keV luminos- 
ity Lo.5-2 keV or higher. Bottom panel: the observed number of 
galaxies in each Lo.5-2 keV bin hosting an ULX of luminosity 
Lo.5-2 keV or higher. 
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39.0 39.5 40.0 40.5 41.0 

Fig. 9. The observed fraction of galaxies with an off-nuclear 
source with a luminosty of Lo.s-2 keV or greater. The red points 
and associated 1 <x confidence region are from our sample, while 
th e dashed line and the l c confidence region have been obtained 
bv lLehmer et al. I ([2006) from the Chandra Deep Fields. 



we could detect an off-nuclear source of 0.5-2 keV luminosity 
Lo.5-2 keV or larger. The bottom panel of the same figure shows 
the observed number of galaxies in each Lo.5-2 keV bin hosting 
an ULX of luminosity Lo.5-2 keV or larger. In order to derive the 
observed fraction of galaxies with an off-nuclear source, we di- 
vided the values of the histogram in the bottom panel by those 
in the top panel of Fig. [8] The result is shown in Fig. [9] The red 
points are the result of our analysis, and the dashed area is the 
lcr confidence regio n computed usin g the prescriptions for small 
numbers statistic bv lGehrelsl (fl986). For compar ison, we report 
in the same figure also the fractions obtained by iLehmer et al. I 
(2006) from the Chandra Deep Fields (dashed line and lcr con- 
fidence region). These fractions should be considered as lower 
limits due to the limited Chandra spatial resolution that does not 
al low to detect off-nucl ear sources with small offsets (see Fig. 10 
of ILehmer et al. Il2006l) . The agreement between our results and 
the CDFs points is reasonably good above log(Lo.5-2 kev) > 40, 
although our point and the associated confidence contours are 
about a factor of two lower th an, but consistent with, those de- 
rived by ILehmer et al. I d2006l) . In the lower luminosity bins, it 
seems that the two measures are discrepant; however, we do not 
consider this difference highly significant, since the measured 
fractions are consistent at the 2cr level. Also, at the faintest fluxes 
the differences between the two X-ray catalogs used is more se- 
vere. For these faint sources the positional uncertainties affecting 
our sample are larger than for the same sources detected in the 
longer Chandra exposures of the CDFs, and therefore we may 
be missing the faintest ULXs in the sample if their error box is 
consistent with the position of the nucleus. We also note that our 
selection criteria for off-nuclear sources reported in Se c. |2] are 
more conservative than the ones used by Lehm eret al. I {2006). 
From Fig.|9]we found that « 0.5% and « 0.2% of the galaxies are 
hosting a ULX with Lo. 5 _ 2 kev > 3x 10 39 and Lo.5-2 kev > 2x 10 40 
erg s~ ! , respectively. 

We now discuss the observed trend of the fraction of ULX 
as a function of th eir X-ray luminosities in the frame of the 
beaming model o f lKingl J2009>. According to this model, ULX 
are stellar mass black holes accreting at a super-Eddington rate 
(m = M/rjLEddC 2 > 1, for a typical radiative efficiency 77 ~ 0.1 
and accretion rate M). Matter accreting at such rates is eas- 
ily blown away close to th e inner edge of the accretion disc 
(Sha kura & Sunvaev"lll973l) : then, the radiative output from the 
resulting flow pattern is of the order of L w LeddU + l n '^ but 
emerges collimated by the central funnel with a beaming^ fac- 
tor b oc mT 2 , so that an external observer who happens to have 
its line of sight within the beaming cone would infer a spherical 
luminosity: Lulx - 10 39 /«7(1 + lmh)/b erg/s ( where mj is the 
black hole mass in units of 7 solar masses; see King! 120091 for 
further details). Thus, neglecting the weak logarithmic depen- 
dence on m, this model directly links the observed luminosity of 
a ULX with its beaming factor b. 

Let us now consider a population of ULX with host galaxy 
space d ensity (as a function o f distance d): n g {d) Mpc~ 3 . The re- 
sults of ILehmer et al. I (|2006) imply an almost linear decline of 
the cumulative number of ULX per galaxy with observed lumi- 
nosity, f bs - fo(iuLx/10 39 ) -1 , where Fq 0.1 is the observed 
fraction of galaxies hosting a ULX with Lulx > 10 39 . The dif- 
ferential fraction <l> bs, i.e. the fraction of galaxies containing a 
ULX with luminosity Lulx per unit logarithmic interval of lu- 
minosity can be derived by simply differentiating the above ex- 
pression, to obtain <J> b s = dN/dLogL\jLX = Lo(Lulx/10 39 )~' ~ 



5 Note that here 'beaming' simply means geometrical collimation, 
and not relativistic beaming. 
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Table 1. Properties of ULXs in C-COSMOS 



xnr 


RA Dec 


Counts^ 




Pos. error^ 


Offset 


Offset 


Offset^ Off-axis£ 




(J2000£ 


(0.5-7 keV) 


(erg s" 1 ) 


(arcsec) 


(arcsec) 


(kpc) 


(arcmin) 



1151 


10:00:10.39 


02:09:23.40 


28 


4O.440.5 


0.5 


1.67 


2.91 


3.42 


2.8 


1388 


10:01:08.46 


02:01:06.05 


17 


40 8^ 


0.6 


3.62 


12.12 


5.73 


2.6 


1870 


10:01:03.76 


02:30:50.22 


9 


39 9^ 

40 6^ 


0.4 


3.12 


4.66 


7.25 


2.6 


2418 


10:00:08.43 


02:14:47.65 


6 


0.3 


1.58 


6.76 


4.53 


1.3 


3441 


09:59:33.78 


01:49:06.92 


5 


40 2 40 5 


0.5 


0.95 


2.26 


1.87 


3.7 


11100 


10:00:58.65 


02:11:39.90 


12 


40 I 403 


0.4 


0.92 


1.85 


2.37 


3.4 


11938 


10:00:43.02 


02:00:32.74 


7 


0.8 


1.39 


4.79 


1.77 


4.4 



" ID of the Chandra source (Elvis et al. 2009) 
b X-ray counts in the [0.5-7] keV band. 
Logarithm of the [0.5-7] keV X-ray luminosity. 
d X-ray positional error. 

e Ratio of the distance between the X-ray centroid and the optical 
centroid over the radius of the Chandra positional error circle. 

t Off-axis angle value in the image where the source is closer to the 
on-axis position. 

g X-ray coordinates of the ULX. 



F(jb/mj, where the last approximate equality has been derived 
neglecting the logarithmic dependence of Lulx on b. 

We now consider the application of such a model to a multi- 
wavelength survey like COSMOS. We define the limiting flux of 
the survey in the X-ray band as /i; m = x 1CT 16 erg s _1 cm -2 , 
so that an object of beaming factor b can be seen out to a dis- 
tance of d{b) = (10 39 m 7 /47r/i iln fe) I/2 313(b/m 7 )- 1/2 f;ll 2 Mpc, 
and express, in full generality, the number density of galaxies as 
a function of distance as n g (d) - n g (d(b)) = n g fl(b/mj) a , where 
n g fl is the number density of possible host galaxies in the survey 
at the maximum distance where an un-beamed source (b — 1) 
can be seen. Such an expression is a very general form appropri- 
ate for power-law luminosity functions in Euclidean Universes 
and is adopted here for the sake of simplicity the exponent a 
depends both on the galaxy luminosity function slope and on the 
survey selection function and can in principle be derived empiri- 
cally for any given survey: typically, for flux-limited surveys, we 
have a > 0, while volume limited ones have a « 0. Given the 
observed cumulative fraction ,F bJZL and the corresponding dif- 
ferential <D b s = Fob/mj one has to search through a space vol- 
ume V ~ l/n g (d)Fo(b/mj) to find a ULX with beaming factor 
b (within a unit logarithmic interval of b). From this expression 
for the volume we derive: 

/ n g0 F \- 1/3 /M" (1+ar)/3 
d(b)= 125 [A — — - (4) 
v ; \ 0.05 0.1/ \m) v ; 

Thus, the minimum beaming factor (corresponding to the max- 
imal luminosity) of a ULX in a survey of area A (in units of 
square degrees), is given by: 

6 Although we have applied a k-correction to the luminosity values 
in equation (1), we resolved to make the calculations in this paragraph 
under the assumption of a Euclidean Universe to simplify the derivation 
of equation (5). 

7 We assume in this calculation that the fraction of galaxies hosting a 
ULX does not change as a function of distance. This is an approxima- 
tion, since the star formation rate varies with redshift and therefore it is 
plausible that the fraction of galaxies hosting a ULX varies too. 



where y = j^rf- The overall efficiency of finding ULX scales as 
InLuLx.max- We recall that Lulx is the spherical luminosity that 
would be inferred by an external observer who happens to have 
its line of sight within the beaming cone. Applying this rough 
estimate with a = 1, 1117 = 1, n ? o = 0.05, Fo = 0.1 to the 
COSMOS survey (/_i 6 2, A = 6.9), we obtain b min - 0.04, 
AjLX,max — 2.5 x 10 40 in reasonable agreement with the present 
data. Interestingly, this also suggests that larger, but shallower, 
surveys could be more efficient in finding ULX (provided a sim- 
ilarly deep sample of host galaxies can be identified): the all sky 
eROSITA survey (/_ 16 = 100, A = 4 x 10 4 ) could find a large 
number of ULX, including microblazars up to b m \ n = 2.6 x 10 -6 , 

LuLX.max - 4.0 X 10 44 . 

5. Conclusions 

We have presented a sample of ultraluminous X-ray sources 
(ULXs) selected from the Chandra survey in the COSMOS area 
(C-COSMOS). From 1761 X-ray sources detected with a maxi- 
mum likelihood threshold of detml=10.8 in at least one detection 
band, we have selected 7 ULX candidates covering the redshift 
range z=0.072-0.283. 

Taking advantage of the excellent ancillary data available in 
the COSMOS field, we have studied the properties of their host 
galaxies. From a detailed morphological analysis of the ACS im- 
ages and rest-frame colours, we found that ULXs are hosted both 
in late and in early type galaxies, with a slight preference for the 
former. 

From the multi-band photometry and from the optical spec- 
tral lines, we have measured stellar masses and star formation 
rates for the host galaxies. Using literature X-ray luminosity 
functions for HMXBs and LMXBs, we have defined probabil- 
ity areas for having detectable off-nuclear sources in the plane 
SFR versus M*. All our ULXs candidates are hosted in galax- 
ies for which we expect a large number of X-ray binaries with 
L > 10 38 erg s -1 . 

The presence of IMBHs (~ 10 2 - 10 5 M Q ) in some of our 
ULXs cannot be excluded with the current data. The best can- 
didates for this new class of accreting black holes are the ULXs 
hosted in early type galaxies (therefore not associated with re- 
cent star formation activity) and with X-ray luminosity above 
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Table 2. Properties of the host galaxies of ULXs in C-COSMOS 



XID^ i mag R P ^ 


?J_ 


Class rf log(MJ 


SFR 


(AB) (arcsec) 




(M ) 


(M /yr) 



1388 
1870 
2418 
3441 
11100 
11938 



15.62 
20.39 
18.40 
19.45 
18.19 
18.48 
18.94 



19.89 


0.094 p 


ETG 


5.20 


0.204 p 


LTG 


9.21 


0.072 s 


LTG 


3.99 


0.283 s 


ETG 


5.50 


0.133 s 


LTG 


6.90 


0.110 s 


LTG 


2.53 


0.221 p 


ETG 



io r 02 

q q+o.2 

-OA 

10.9$ 
10.5$ 

10 4 + ° 1 

11 
-o.i 



3.0+ 08 



10.9 + 



0.3_ 
l.l+o 

i.7 + ; 

0.04+! 



" ID of the Chandra source (Elvis et al. 2009) 

* Petrosian radius of the host galaxy Petr osianl i 19761) 
c Redshift of the host galaxy: "s" for spectroscopic and "p" for pho- 
tometric redshifts. 

d Morphological classification of the ULX host galaxy: early type 
galaxy (ETG) or late type galaxy (LTG). 

e The photometric coverage is limited to few bands and we cannot 
constrain its M, or SFR 



10 41 ergs" 1 that can be difficult to explain with high-mass stellar 
black holes. The objects that satisfy these criteria from our sam- 
ple are XID= 2418 and 11938. Longer X-ray exposures could 
give us more insights on the real nature of these sources from 
a detailed study of the X-ray spectrum. Similarly, we cannot set 
constraints on the recoiling black-hole nature of our sources with 
the curr ent data, but it is worth m entioning that recent predic- 
tions by Volon teri & Madaul d2008l) expect at most one of such 
objects in the C-COSMOS survey, assuming the most favorable 
scenario (spinning black holes, no bulge in the host galaxy, long 
active phase). 

Finally, we have derived the fraction of galaxies hosting 
a ULX as a function of the X-ray luminosity. We found that 
ss 0.5% and a 0.2% of the galaxies are hosting a ULX with 
L .5- 2 fev > 3 x 10 39 and L . 5 -2fev > 2 x 10 40 erg s" 1 . This 
is in reasonably good agreement wi th the observed fract ion de- 
rived in the Chandra Deep Fields by Lehme r~et al. I J2006) above 
log(Lo.5-2fcv > 40 erg s . A possible discrepancy in the lower 
luminosity bins can be likely attributed to the differences in the 
limiting fluxes of the two catalogs and, therefore, to the different 
positional uncertainties affecting faint X-ray sources. 
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